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I ' Abstract 

^! In view of the imperfection of the Navier-Stokes equations in wave de- 

fj I scription, a simple new model is developed here. Its expression differs from 

£^ ■ the common trigonometric form in that: (1) it is in a coordinate form syn- 
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thesized by two piecewise quadratic-curve functions; (2) it has more distinct 
physical meaning; (3) it accords well with the actual one which has a sharp 



>■ , crest and a flat trough; (4) it yields new relations for the wave parameters, 

m 

^ \ which indicate the increasing of amplitude implies the decreasing of frequency 



with a quicker speed. These understandings may promote the reconstruction 



(^ ! of ocean wave theory. 
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1. Introduction 

The study of ocean wave is the oldest branch of hydrodynamics. If it is 
timed from the apphcation of potential theory the duration should be no less 
than 200 years and now the classical wave theory is developed, such as in 
[1-4]. Certainly, the default model is the Navier-Stokes equations, after all, 
the ocean wave is involved in fluid pressure and movement velocities. But all 
in all, the Navier-Stokes equations are very complex nonlinear ones and the 
classical results for ocean wave are merely some kind of linear approximations. 

On the assumption that the amplitude A of a wave on the water surface 
is infinitely small relative to the wave-length A, that is, A/\ <^ 1, there is a 
classical linear wave theory. At this time, the surface traveling wave can be 
expressed in the trigonometric formt^"^': 

C,{x,t) = Asm{kx — ut), (1) 

where k and u are the wave number and frequency. By the way, these 
two parameters are the inverse expressions for the wave-length A and wave- 
period T with k = 27r/A and u = 2tt/T. For the deep-water case, there is a 
dispersion relation: 

u' = gk (2) 

and a wave-speed formula: 

^ ^ fg 9 /ON 

^=T = I = V"-' ^'^ 

here g is the gravity acceleration, which means a wave with lower frequency 
or wave-number should travel faster. 



In case A/\ (= Ak/2Ti) is not infinitely small, that is, the wave steepness 
Ak takes a finite small value (< 1), there is a finite-amplitude Stokes wave 
theory'^'^1. At this time, the wave surface and dispersion relation for the 
deep-water case become 

{(x, t) = Acose + Q.hkA^ cos 26 + Q^hk'^A^ cos 3^ + ■ ■ ■ , (4) 
''-gk[l + {Akf + l.2h{AkY + •■■], (5) 
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where 9 = kx — ut. Moreover, the wave speed can be also deduced from 
c = oo/k. By the way, for a too big steepness Ak which accords with a too 
sharp crest the Stokes wave would break down^^l. 

The above results are very classical and they take fundamental parts 
in the whole gravity wave theories. As the name implies, a wave is called 
"gravity wave" if its main restoring force is gravity. We note that, besides 
the gravity, the maintaining of up-down fiuctuation also needs an upward 
buoyancy. Do the frequency u, wave-number k and wave-speed c depend on 
the amplitude A? This is an open problem for the above classical wave 
theories, though the dispersion relation for Stokes wave includes the effect 
of amplitude. After all, the wave-number k is also an unknown parameter, 
the dependence of frequency u on amplitude A is still a suspense, not to say 
the wave-speed c. As the problem concerned there is not much hope for the 
Navier-Stokes equations. The reason is that, the solving of these equations 
involves in determining the upper surface boundary condition which is just 
the wave to look for^. In fact, even the trigonometric wave forms themselves 
are some kind of conjectures. So it needs an innovation on the wave model. 



Notice that the ocean wave is a synthesis of transversal wave and longi- 
tudinal wave, we decompose the motion into the vertical part and horizontal 
part and model them separately in the next two sections. The motion synthe- 
sizing process and the traveling-wave modeling process are given in Section 
4 and 5 respectively. These modeling processes result in explicit relations for 
the wave parameters in Section 6 which differ much from that of the classical 
wave theories. To reflect the characteristics of the derived model we have 
also made some numerical simulations for it in Section 7. In addition, the 
summary together with some related discussions are given in the last section. 

2. Modeling the Particle's Vertical Motion 

We assume that the particle at the water surface always maintains at 
it during the reciprocating motion process. In fact, this is a fundamental 
assumption for wave study, such as in [3]. To simplify the problem only 
the deep-water case is considered here. We note that the small-amplitude 
assumption is abandoned and the effect of surface tension is neglected. Notice 
that the particle's horizontal motion has no effect on the vertical one we 
simplify the problem to an one-dimensional one. 

Take the undisturbed water surface as the zero-level {z = 0) with positive 
displacement z{t) upward. The corresponding vertical velocity is denoted by 
w{t). To be simple, the start time t = is chosen from the undisturbed 
state, that is, z{0) = and the particle's upward motion can be understood 
as from a sudden impulse. In the following we model the reciprocating motion 
process by 4 stages. 



2.1. Stage 1: from the Zero-Level to the Crest 

Due to the sudden impulse, the surface particle gains an initial upward 
velocity w(0) = wq. In the following we describe the motion by force analysis. 
Denote m and a^ the mass and vertical acceleration of the particle separately. 
It follows from the Newton's second law F = ma that 

maz = -mg + P«, - Pa, (6) 

where Pw and Pa stands for the upward pressure from the water below and 
the downward pressure from the air. In case there is no wind or the wind is 
not so strong, we can neglect the air pressure effect and take Pa = 0. As for 
the upward pressure P^, though it plays a part during the rising process, it 
should be smaller than the gravity mg, after all, this is a decelerating process 
and the velocity should be at the crest. Roughly, we can neglect the effect 
of Pyj and take the gravity m,g as the main restoring force. This simplicity 
accords well with the concept of "gravity wave". Under the above analysis, 
we can approximate the vertical acceleration by the gravity acceleration, that 
is, a^ ^ -g. 

At this time, the upward velocity and displacement satisfy 

wit) = Wo - gt, (7) 

zit) = Wot - ]^gt\ (8) 

Particularly, when the surface particle attains the crest with w{ti) = and 
z{ti) = A, these equations yield 

In A 

(9) 




where A is actually the wave amplitude. 

2.2. Stage 2: from the Crest to the Zero-Level 

For the falling case, the eqn. (6) and its approximation a^ ~ —g also hold. 
At this time, the motion can be seen as in a free-falling manner with reset 
initial time t = at the crest. So the downward velocity and displacement 
satisfy w(t) = —gt and z{t) = A — gt'^/2. Particularly, when it attains the 
zero-level with z{t2) = its time is 

t2 = J-. (10) 

V ^ 

2.3. Stage 3: from the Zero-Level to the Trough 

It follows from Stage 2 that the particle possesses an initial velocity 



w{0) = —gt2 = —\/2Ag at the zero-level with a reset time t = 0. At 
this stage the falling goes on with a decreasing velocity due to the existence 
of upward buoyancy. To do the force analysis we choose a tiny water col- 
umn as the object. Its volume is defined in the form V(t) = h{t)S, where 
h{t) = —z(t) is the distance between the zero-level and the water surface 
and S is the horizontal area of the water column. For this case, the buoy- 
ancy Fi){t) is determined by the volume of water column together with that 
of the emptied air space in the same quantity. So Fh{t) = 2pgV(t), where 
p stands for the water density. Notice that the mass and gravity equal to 
m(t) = pV{t) and G{t) = pgV{t) separately, it follows from Newton's second 
law that Fi,{t) — G{t) ^ m{t)az{t) which results in az{t) ~ g. This indicates 
that the particle at the surface possesses an upward acceleration g. So the 



falling velocity is w{t) = w{0)+gt = —gt2 + gt and when it attains the trough 
with w{t3) = the cost time should be ta = ^2 = a/2A/(7. 

2.4- Stage 4- from the Trough to the Zero-Level 

This stage is the inverse process of Stage 3. At this time, the particle 
accelerates upward from the trough. When it reaches the zero-level, it follows 
from = z{t4) = —A + gt\/2 that t4 = \J2A/ g. Moreover, it is easy to see 
that w{ti^ = gt4 = gti = Wq which implies the beginning of anther up-down 
period. 

By combining the previous four stages we get the particle's up-down re- 
ciprocating motion period: 

l2A 

T = ti+t2 + t3 + t4=4W — . (11) 

To arrange the four stages with a uniform time, then the vertical displacement 
can be rewritten in a piecewise quadratic-curve form: 

f A-g(t- T/Af /2, < t < T/2, 
z{t) = { ''' ' - - I ^ (12) 

\ g{t- ?,T/Af /2 - A, T/2<t< T. 

3. Modeling the Particle's Horizontal Motion 

Besides the vertical motion, there is also a horizontal one for the surface 
particle. This motion is mainly caused by the horizontal pressure-gradient 
force at the slant water surface. 

For a hydrostatic case, the water pressure P increases linearly along with 
the water-layer thickness s, that is, P = pgs. Assume that the "hydrostatic 
approximation" holds for the wave fluctuation case, then P ^ pg[^{x, t) — z]. 



here $,{x,t) and z are the surface elevation and water depth below the zero- 
level (with positive sign upward). Particular, for the surface particle it ac- 
cords with the case z = 0. Under this assumption, the horizontal acceleration 
can be approximated by the Euler equation'^'^': 

a, = ^ = -1^ = -g^ (13) 

^ dt p dx dx 

To simply the problem, we take the wave slope d^/dx as a parameter and 
denote 5 the average of its absolute value over a wave-length, that is, 

(9^(a;,to) 
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where to is a fixed moment. With this simplicity, the horizontal acceleration 
can be further approximated by ax ~ 5g for the case d^/dx < and ax ~ —5g 
for the case d^/dx > 0. In the following we describe the particle's horizontal 
motion according to the 4 stages adopted by the vertical modeling. 
3.1. Stage 1: from the Zero-Level to the Crest 

We choose the start time t = as the moment that the surface particle is 
just at the zero- level. Alone with the rising movement, it accelerates horizon- 
tally due to the existence of positive pressure-gradient force with dC,/dx < 0, 
after all, the crest is on the left side of it at this time. So the horizontal ve- 
locity satisfies u{t) = uo-\-6gt, where uq is the initial velocity which is usually 
caused by the current. Particular, for the no-current case it reads u{t) = 6gt. 
Correspondingly, the horizontal displacement relative to the location of the 
crest is 

x{t) = -6A+-5gt'^, (15) 

since the crest occurs at t = T/4 with vertical displacement A = (?(T/4)^/2. 



3.2. Stage 2: from the Crest to the Zero-Level 

When the particle attains the crest its vertical velocity decreases to and 
its horizontal velocity increases to the maximum value Uc = 6gT/4. So the 
next stage is a falling process. At this time, the crest is on the right side of it 
with d^/dx > 0. So u{t) =u^- 6g{t - T/4) with u{T/2) = 0. Similarly, the 
corresponding horizontal displacement relative to the location of the crest is 

x{t)=u,{t-T/4)-^6g{t-T/4f, te[T/4,T/2]. (16) 

3.3. Stage 3: from the Zero-Level to the Trough 

At this stage, the falling goes on with a negative horizontal acceleration. 
So u{t) = 0—6g(t—T/2). Relative to the location of the crest its displacement 
becomes 

x{t) = 6A-^6g{t-T/2)\ t e [T/2,3T/4]. (17) 

3.4. Stage 4'- from the Trough to the Zero-Level 

When the particle attains the trough its vertical velocity decreases to 
and its horizontal velocity decreases to the minimum value Ut = —6gT/4. 
As it increases, the horizontal acceleration turns to positive which results in 
u{t) =ut + 6g{t- 3T/4) and 

x{t) =ut{t- 3T/4) + ^6g {t - 3T/4)' , t G [3T/4, T]. (18) 

Moreover, it is easy to see that x{T) = —6A and u{T) = which imply the 
beginning of anther wave period. 



To synthesize the above 4 stages with some modulations, it yields a model 
for the surface-particle's horizontal motion: 

{-5 A + 5gt^/2, < t < T/4, 
6A~6g{t-T/2f/2, T/4<t<3T/4, (19) 

-5A + 5g{t - T)V2, 3T/4 <t<T. 

Here the default equilibrium point is xq = 0. Similar to the vertical model, 
this is also a piecewise quadratic-curve form. 

4. Synthesis of the Particle's Motion 

According to the previous modeling processes, the trace of the surface- 
particle's motion can be seen as the location shifting of {x{t),z{t)). Here 
x{t) and z{t) are the horizontal and vertical coordinates given in eqn.(19) 
and eqn.(12) separately. It follows from Figure 1 that the models of piecewise 
quadratic-curve forms for the horizontal and vertical motions are very close 
to the trigonometric forms adopted by the classical linear wave theory with 
some modifications. This, at least, indicates that the derived models are 
acceptable. Their advantages lie in distinct physical meaning, flexible wave 
slope and explicit dependence of frequency on amplitude. In addition, as 
shown in Figure 2, the trace of the particle's motion is very close to the 
elliptic curve x'^/{6Ay + z'^/A'^ = 1. It is interesting to find that the lower 
wave slope 6 implies the thinner ellipse. 

5. Modeling the Traveling Wave 

Since the wave fiuctuation can be seen as a coherent movement of many 
surface particles, it is convenient to express the wave surface in a coordi- 

10 
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Figure 1: The variations of the horizontal and vertical displacements for the surface particle 
with A — S, 



6 = 1/3 and w = ir^/g/SA. 
Asinujt and —5Acos{ujt) are also drawn. 



For comparison, the trigonometric functions 
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Figure 2: The trace of the surface-particle's motion with A = 3, S = 1/3 ( — ). For 
comparison the elliptic curve x'^ + z'^ /9 = 1 (- -) is also drawn. Here the arrows depict the 
movement velocities at the key points. 
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nate form {x{xQ,t),z{xQ,t)). Here Xq is the common horizontal location and 
x{xo,t), z{xQ,t) are the horizontal and vertical displacements of the surface 
particle relative to the equilibrium xq. To accords with the using habit, some- 
times we still use the nominal denotation z = C,{x,t) as in the previous to 
stand for the wave surface, though x and z have no explicit relation. At first, 
we describe the instantaneous state of wave surface at t = in the following 
way. 

To incorporate the horizontal location xq into eqn.(19) by a substitution 
t = kxo/u together with a translational term, it leads to: 

-6A + 6g{kxof/2uj'^ + xo, < A;a;o < 7r/2, 
x{xo,0)= { ^A-^9ikxo-Trf /^u^ + xo, n/2<kxo<3n/2, (20) 

-5 A + 6g{kxo - 27r)V2a;2 ^ ^^^ 3^/2 < kxo < 2-k. 

Here the transform T = 2ti /oo is used. For the vertical one, since our modeling 
is started from the moment that the surface particle is at the zero-level with 
an rising trend, the wave surface should begin with a trough rather than a 
crest. So it follows from eqn.(12) that 

^ ( -A + g {kxo - 71 /2f /2u^ 0<kxo<7r, 
z{xo,0) = < „ (21) 

[ A- g {kxo - 37r/2)^ /2u^, n < kxo < 2tt. 

Now that the initial wave surface is known, the traveling wave form only 
requires a substitution of kxQ with 6 = kxQ —ut — 2n7r (n = 0, 1, 2, ■ ■ ■). The 
final traveling wave of coordinate form is (a;(xo, t), z(xo, t)) with 
-5 A + bgQ'^l2u? + xq, < ^ < 7r/2, 
a;(a;o, t) = \ ^A-bg{e- ixf /2uj^ + xo, 7^/2 <6< 37r/2, (22) 

-5A + 5g{e-2Tif/2uj'^ + xo, 37r/2 < ^ < 27r. 

12 



z{xo,t) 



(23) 



-A + g{e-Tr/2f /2uj^, < 9 < tt, 
A-g{e-3n/2f /2uj'^, n<9<2n. 
These are piecewise quadratic-curve forms which differ much from the trigono- 
metric form given in eqn.(l). 

6. Formulas for Wave Parameters 

It follows from eqn.(22) and (23) that the wave parameters associated 
with the derived traveling wave are A, u, k and 6. To derive the relations 
for them we take the wave amplitude A and the averaged wave slope 6 as fun- 
damental ones. Firstly, eqn.(ll) yields an explicit dependence of frequency 
on amplitude: 



u 
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(24) 



Notice that the wave slope keeps its sign on each A/4 with almost the same 
averaged value, set to = then it follows from eqn.(14) that 
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which yields 



Tl5 

~ 2A' 



Further deduction leads a formula to the wave speed: 



A oo 



25 ' 



(25) 



(26) 



which indicates that the wave with a big amplitude or a small slope should 
possess a quick speed. To combine u and k with A it leads to a dispersion 
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relation as follows: 



.2 ^ 



^' = -^^9k. (27) 



This indicates that the one for linear wave accords with the case 5 = 7r/4 ^ 
0.78 which is a too sharp slope to meet the small-amplitude hypothesis. In 
another word, the dispersion relation cu^ = gk for ocean wave is doubtful. Due 
to the nonlinear approximation processing on the Navier-Stokes equations, 
the one for Stokes wave should be more reliable. By substitution of wave 
steepness with Ak = it6/2, the dispersion relation in eqn.(5) for the Stokes 
wave can be rewritten as w^ = aogk with ao = l + (7r(5/2)^ + 1.25(7r5/2)'^ + - ■ •. 
So for the actual wave, its dispersion relation w^ = agk should possess a 
factor a > 1 and our one is valid for 6 < 7r/4. As for more exact upper bound 
of 5 it would involve in wave-breaking arguments^^'^^ and further research is 
expected. 

7. Numerical Simulations 

To reflect the characteristics of the derived model we make some numerical 
simulations here. Chose A = 3 ra and 5 = 0.6, then to = Tx^fgjKA ^ 2.00 
rad/s and k = 7t6/2A ^ 0.31 rad/m. It follows from Figure 3 that the wave 
of piecewise quadratic-curve form accords well with the actual one which 
has a sharp crest and a fiat trough. For comparison, the linear wave of 
trigonometric form —A sin {kx — ut) and the third-order Stokes wave at t = 
are also drawn. The differences are very clear. By the way, here only the 
case t = is considered and there is no reflection on the frequency diversity 
which can be tested with further simulations. In addition to the surface state 
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at t = s, the one at t = 0.7T = {).7y^^2A/g ^ 2.12 s is also drawn for our 
model. It indicates that the traveling characteristic accords well with the 
common sense. 




o 



20 30 40 

Distance (m) 

Figure 3: Numerical simulations for different wave models, a: the linear wave of trigono- 
metric form —Asm{kx — cut) at i = 0; b: the third-order Stokes wave at t = 0; c: the 
wave of piecewise quadratic-curve form at i = (— ) and t = 0.7T (- -). Here k — ttS/2A 
with A = 3 and 6 = 0.6. 



8. Summaries and Discussions 

By modeling approach we have deduced a new model for the ocean wave 
which can be seen as an innovation on the classical ones. It differs from the 
common trigonometric solution of Navier-Stokes equations in that: (1) It is 
in a coordinate form synthesized by two piecewise quadratic-curve functions; 
(2) It has more distinct physical meaning; (3) It accords well with the actual 
one which has a sharp crest and a flat trough; (4) It yields new relations for 
wave parameters: 



I 9 , TiS 

uj = 7r\ — 7, k = — 7, 
8A' 2A' 
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V2^ 
26 ' 



which indicates the crack of the open problem. These formulas are very 
useful. Here we only give two examples for the first one. 

If the wave is under a developing process, the first relation can give a 
well guidance to its variation. That is to say, along with the increasing of 
amplitude A, its period increases in the clear way with T = sJ'ilAj g k, 
\.%\/A. It follows from Figure 4 that if the wave amplitude is increased from 
4 m to 5 m then the corresponding period should increase from 3.6 s to 4 s. 




3 4 5 6 
Amplitude (m) 



Figure 4: The variation of wave period T along with the increasing of amphtudc A. 



This formula can be also used in a reverse way. Notice that the amplitude 

is A = g-n"^ j%iJ^ ^ we can express the wave energy by 



^ = 1.4^ = ^.!!. 

2 128^4 ^4 



Therefore, the wave with a low frequency should possess a high energy and 
the dependence is up to the fourth power. This understanding is helpful for 
improving the wave spectrum theory'^l 

In addition, in view of the similarity in Figure 1, the ocean wave model of 
piecewise quadratic-curve forms (22) and (23) can be approximated by the 
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simple fully-smooth trigonometric forms: 

x{xo, t) = —6A cos{kxQ — ut) + xq, (28) 

z{xo, t) = —A sin(A;a;o — ut), (29) 

which takes the linear wave z = —Asin{kx — ut) as a limiting case with 
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